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ABSTRACT 

The isotropy of the Lyaforest in real-space uniquely provides a measurement of cosmic geometry at z > 2. 
The angular diameter distance for which the correlation function along the line of sight and in the transverse 
direction agree corresponds to the correct cosmological model. However, the Lya forest is observed in redshift- 
space where distortions due to Hubble expansion, bulk flows, and thermal broadening introduce anisotropy. 
Similarly, a spectrograph's line spread function affects the autocorrelation and cross-correlation differently. In 
this the second paper of a series on using the Lya forest observed in pairs of QSOs for a new application of the 
Alcock-Paczyriski (AP) test, these anisotropics and related sources of potential systematic error are investigated 
with cosmological hydrodynamic simulations. Three prescriptions for galactic outflow were compared and 
found to have only a marginal effect on the Lya flux correlation (which changed by at most 7% with use of the 
currently favored variable-momentum wind model vs. no winds at all). An approximate solution for obtaining 
the zero-lag cross-correlation corresponding to arbitrary spectral resolution directly from the zero-lag cross- 
correlation computed at full-resolution (good to within 2% at the scales of interest) is presented. Uncertainty 
in the observationally determined mean flux decrement of the Lya forest was found to be the dominant source 
of systematic error; however, this is reduced significantly when considering correlation ratios. We describe 
a simple scheme for implementing our results, while mitigating systematic errors, in the context of a future 
application of the AP test. 

Subject headings: cosmology: miscellaneous — intergalactic medium — methods: numerical — quasars: 
absorption lines 
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1. INTRODUCTION 

Significant observational and theoretical advances in re- 
cent decades have made the Lya forest a powerful and 
unique cosmological tool for studying the hi gh-redshift uni- 
verse. Originally named dWeymann et al.l[l981l) for the dense 
pattern of seemingly discre te Lya abso rption lines seen in 
high-redshift QSO spectra (Lvndsl [l971l) . the absorption is 
now understood to trace a continuous distribution of non- 
uniform neutral hydroge n gas that in turn maps the un- 
derlying dark matter (see lRauchl [l998 for a review). The 
competing processes of recombination and photoionization 
lead to a tight relationship between the density of the gas 
and the neutral fraction, giving rise to a relatively straight- 
forward link between Lya absorption and the large scale 
structure of the universe. Cosmological simulations em- 
ploying this prescription have had remarkable success re- 
producing detailed properties of the Lya fo rest provided by 
high-resolution ground-based QSO spec tra (Cen et al.] ll994t 
IZhang et allll995t iHemarast et alJll996t iTheuns et al.lll998l) 
and lo w-z HST observations dPetitiean et al.lll995UDave et alj 
fl999h . paving the way for the Lya forest to be reliably used 
f or cosmological in vest igation. 

iHuiet al.ldl999l) and lMcDonald & Miralda-Escudel (1 19991) 
first suggested using autocorrelation and cross-correlation 
measurements in the Lya for est for a new applicat i on of the 
Alcock-Paczyhski (AP) test (lAlcock & Paczynskil 1 19791) . a 
purely geometric method for measuring cosmological param- 
eters that is primarily sensitive to J7a at z > 1 . The essence 
of this cosmological test is that spherical objects observed at 
high redshift will only appear to be equal in their radial and 
transverse extent if the correct angular diameter distance is 
used to determine the latter. More generally, the correlation 
function of an isotropic medium, such as the Lya forest, mea- 



sured as a function of separation along the line of sight (the 
autocorrelation £||) and in the transverse direction (the cross- 
correlation will agree only if the correct cosmology is 
assumed. 

Spectroscopy of the Lya forest in any single QSO spec- 
trum yields the complete autocorrelation, albeit with signif- 
icant variance from one line of sight to another. The cross- 
correlation, on the other hand, must be pieced together from 
pairs of QSOs with different transverse separations. Until re- 
cently, only approximately a dozen pairs with similar redshifts 
(so that their Lya forests overlap) and separations of a few 
arcminutes or less (the correlation signa l diminishes rapidly 
beyond this point) were known (see, e.g., Rollin de~et al.ll2003l 
and references ther ein). The 2dF QSO Redshift Survey (2QZ; 
ICroom et al.ll2004) s ignificantly increa sed this number, and 
with motivation from Mc Donald! (12003). moderate resolution 
spectra (FWHM ~ 2.5 A) with modest signal-to-noise ratios 
(S/N > 10 per pixel) hav e been obtained for m ore than 50 
of these p airs at the VLT dCoppolani et al.ll2006h . MMT, and 
Magellan (Mar ble et al.ll2.Q08L hereafter referred to as Paper I) 
observatories. 

While conceptually simple, the Lya forest variant of the AP 
test is not as straightforward as measuring the autocorrelation 
and cross-correlation from pairs of QSOs and determining the 
angular diameter distance which satisfies the presumption of 
isotropy. Rather, two additional sources of anisotropy must 
be accounted for. First, the line spread function (LSF) of the 
spectrograph smooths QSO spectra along the line of sight, af- 
fecting the autocorrelation and cross-correlation differently. 
Second, nonzero velocities caused by the expansion of the 
universe, gravitational collapse, and thermal broadening make 
the correlati on function in redshift-space (z-space) a nisotropic 
dKaiserlll987l) . Fortunately, the theoretical work of iHui et ail 
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1999). lMcDonald & Miralda-Escudeld 19991) . and lMcDonaldl paper. 



(2003) found that these redshift-space distortions can be dis- 



entangled from the desired cosmological signature. 

The aim of this paper is to investigate these non- 
cosmological anisotropics in the Lya flux correlation func- 
tion in a manner which is directly applicable to observations 
of QSO pairs that are suitable for a new application of the AP 
test. To this end we have employed a variety of cosmological 
hydrodynamic simulations to model the autocorrelation and 
cross-correlation of the Lya forest in both redshift-space and 
real-space. These simulations are described in § [2] as well 
as our procedure for ex tract ing mock Lya absorption spectra 
from them. In §[3]and § !3.1l we introduce the correlation func- 
tion and discuss how to mitigate relevant size and mass res- 
olution limitations of the current generation of simulations. 
The effect of arbitrary spectral resolution on the correlation 
function is the subject of § 13.21 Additional potential sources 
of systematic error, both computational and observational, are 
addressed in § H] The implications of our results for the AP 
test are the topic of § [5] Finally, in § [6] we summarize this 
work and its findings. 

2. SIMULATION DATA 

2.1. Simulations 

This body of work draws from a suite of eight cosmologi- 
cal simulations which primarily differ in their size, mass res- 
olution, and prescription for galactic outflow (Table [TJ. To- 
gether, wl6n256vzw (abbreviated as wvzw) and g6 mitigate 
the effects of limitations in volume and mass resolution as 
discussed in § 13.11 Differing wind models (described in § !4.51 l 
are investigated with the wl6n256cw and wl6n256nw simu- 
lations (abbreviated as wcw and wnw respectively), which are 
otherwise identical to wvzw. Similarly, the ql-q4 simulations 
differ only by their number of particles, N p , and are used in 



§ 13.11 to test for convergence as a function of mass resolu- 
tion. All of these simulations have been the subject of previ- 
ous study; therefore, we address only the relevant details and 
direct interested readers to the references provided for addi- 
tional discussion. The common genesis of these simulations 
is described below, while their differences are contrasted in 
Table[TJand the sections referenced above. 



The jV-body + hydrodynamic code G ADGET (Springel et al 



2001b. with modifications described in Springel & Hernquist 
2003), was used to create ql, q2, q3, q4, and g6 



Springel & Hernquistl I2003L but also see iFinlator et al.1 
(2006) regardi ng g6 ), while wvzw, wcw, and 



Oppenheimer & Dave 2006) were run with a similarly mod 



ified version of its successor Gadget-2 (Spri ngell 120051) . 
This code computes gravitational forces via a tree particle- 
mesh solver and hydrodynamical forces with an entropy- 
conservative formulation of smoothed particle hydrodynam- 
ics (SPH). Modelling of additional physical processes in- 
cludes prescriptions for star formation and supernova feed- 
back within evolving galaxies, which impact the intergalac- 
tic medium (IGM) via outflow from galactic winds. A 
spatially uniform photoionization background is included, 
with the spectra l shap e an d redshift evolu t ion gi ven by 
lHaardt & Madaul (1 19961) and lHaardt & Madaul (1200 ll) for the 
Gadget and Gadget-2 runs respectively. Radiative heat- 
ing and cooling is calculated assuming photoionization equi- 
librium and optically thin gas. All of the simulations were 
run as cubic volumes with periodic boundary conditions and 
the same cosmological parameters (f2„, = 0.3, J7a = 0.7, fi/, = 
0.04, c% = 0.9, and h = 0.7), which we assume throughout this 



2.2. Line of Sight Selection 

"Observed" lines of sight through the simulation box (par- 
allel to the three principal axes) were grouped into sets, each 
of which probed the desired separation range (0-5 arcmin- 
utes). N s sets were randomly distributed across each of the 
three mutually orthogonal faces of the box in order to rep- 
resentatively sample the diversity of structure present in the 
simulation volume. The value of N s (Table [U, which roughly 
scales inversely with the box length of the simulation cube 
for comparable total path length, yields oversampled structure 
(correlated measurements) in some cases. However, there are 
sufficient independent lines of sight through each simulation 
to achieve negligible uncertainty in the Lya forest flux corre- 
lation mean despite considerable variance. 

The lines of sight in a given set were arranged in the fol- 
lowing manner. A position on the face of the simulation box 
was randomly selected, from which an imaginary 300" long 
line was extended at a random angle within the same plane. 
If the line happened to intersect an edge of the simulation 
face, it was continued on the opposite side, per the wrapped 
boundary conditions. The two opposing ends of the line and 
11 intermediate positions (2, 3, 8, 12, 15, 25, 45, 80, 150, 
210, and 250 arcseconds from the origin) defined starting co- 
ordinates for that set. The intervals between these 13 lines 
of sight, determined via a Monte Carlo approach designed 
to maximize sampling of angular separation (particularly at 
small separations where the correlation function evolves more 
rapidly) with a minimal number of spectra, yield 73 pairings 
with unique separations. 

2.3. Lya Flux Spectra 

Two Lya transmitted flux spectra, one corresponding to 
redshift-space and another to real-space, were computed as 
a function of position, x, along each line of sight using a 
modified ve rsion of the progra m specexbin (originally part 
of tipsy; iDave et aTll 19991) . First, the physical properties 
of the gas (density, temperature, and velocity) were calcu- 
lated at ~ 20 h' 1 comoving kpc intervals (Av ~ 2.1km s" 1 
or AA ~ 0.029 A at z = 2.4). For the real-space spectra, the 
velocities (resulting from Hubble expansion across the length 
of the box, bulk flows, and thermal broadening) were reset 
to zero. Then, the corresponding H I opacities, r, were de- 
termined, using ioniz ation fraction look up tables generated 
with Cloudy v9 6 dFerland et al.l 1 1 998b . and converted to 
Lya transmitted flux, 



fix)- 



-T(X) 



(1) 



In an additional intermediate step, the extracted opacities 
were multiplied by a single scaling factor in order to match the 
mean transmitt ed flux, (/), in red shif t-space to the observe d 
value of either iPress et al.l d 1993b or iKirkman etail d2005l) . 
For the moderately overdense regions characteristic of the 
Lya forest (p/p < 10), this has the same effect as chang- 
ing the amplitude of the photoionizing background (which 
determines (/)) when running the simulation and/or when 
later computing th e ionization fraction to determine opacity 
dCroft et alj|1998l) . The appropriate scaling was determined 
for each simulation and redshift via an iterative process. The 
raw opacity values from all redshift-space extractions were 
multiplied by a single scaling factor (originally one), con- 
verted to fluxes, and averaged. The scaling factor was then 
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adjusted to be higher or lower as needed, and these steps were 
repeated with incrementally smaller adjustments until the re- 
sulting mean flux agreed with the desired value. This con- 
vergence was considered complete when the difference was 
less than or equal to the formal error in the calculated mean. 
Given the large number of extracted opacities, this typically 
corresponded to a relative difference of less than 0.01 %. 

Figure Q] shows the resulting spectrum for a single line of 
sight through the wvzw simulation box at different redshifts. 
The corresponding set of lines of sight (for z = 3) is shown in 
Figure [2] illustrating the decreasing correlation (in the form 
of visual similarity) at increasingly larger transverse separa- 
tions. A comparison of the same line of sight in real-space 
and redshift-space is provided in Figure [3] where the redis- 
tribution of opacity due to redshift-space distortions is subtle, 
but evident. In addition, narrower absorption features can be 
seen relative to a different line of sight through the larger g6 
simulation, due to the poorer mass resolution of the latter. 

3. THE CORRELATION FUNCTION 

From the ensemble of lines of sight through each simula- 
tion, we know the transmitted Lya flux as a function of veloc- 
ity, 



1 "f" Zsim 



(2) 



in the radial (vn) and transverse (vj_) directions for En = N s 
different realizations (i.e., sets of lines of sight). Here H(z S i m ) 
is the value of the Hubble parameter at the fixed redshift of 
the simulation, and the denominator accounts for x being in 
comoving coordinates. For notational convenience, we define 
S to be the relative difference between / and the global mean, 



<5« (vi|,v_l) 
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The relation between transmitted flux separated along the line 
of sight or in the transverse direction by a velocity difference 
Av is given by the autocorrelation, 

£n (Av) =jrJ2jrJ2]zJ2 6 » ( v \\p v ^) s " ( v \\j + Av ' v ±) • 

' S «=1 i=l II /=1 

(4) 

and zero-lag cross-correlation, 



1 



A', 



^ Av > = ^ Th ^ 6 " ( v ll n S " ( v \\j > V U + Av ) > (5) 

s «=1 II j=l 

respectively. Note that /„ (and therefore S„) is periodic due 
to the wrapped boundary conditions of the simulation box. In 
real-space, the correlation of the Lya forest is isotropic, and 
£|| = Figure|4]confirms this basic result for the hybrid cor- 
relation measurements (discussed in § 13. Il l at z = 3 and shows 
the anisotropy introduced in redshift-space. 

3.1. Box Length vs. Mass Resolution 

The Lya forest is believed to have formed, via gravitational 
collapse, from perturbations in the initial density field. In or- 
der to reliably model the correlation function of the Lya for- 
est, simulations must evolve a sufficiently large volume with 
adequate mass resolution. A simulation box length that is too 



small, or a gas particle mass that is too large, excludes rele- 
vant perturbations on large and small scales respectively. In 
the moderately overdense regime of the Lya forest, growth is 
sufficiently non-linear that perturbations of different sizes be- 
come coupled, and the correlation function is affected even at 
scales not excluded. In addition, aliasing due to the periodic 
boundary conditions of the simulations is extended from half 
the box length (L) to smaller scales (we limit our analysis to 
separations less than L/3-L/4). Since simulations which can 
satisfy both of these competing demands are not yet available, 
we mitigate these effects by forming hybrid correlation curves 
from two different simulations which meet the requirements 
independently. 

The wvzw simulation has 256 3 gas particles within an 
L = 16 h~ l Mpc box, yielding a gas particle mass of m gas = 
2.71 x 10 6 h Mq. In order to verify that this mass reso- 
lution is sufficient for our purposes, we used the ql, q2, q3, 
and q4 simulations to test for convergence (note that the result 
may be simulation code dependent). The ^-series are identi- 
cal except for particle number, with gas particle masses which 
decrease with increasing series number (42.4, 12.5, 3.72, and 
1.10 in units of 10 6 h~ l M ). As a consequence of their small 
box length (L = 10 h~ l Mpc), the correlation is artificially de- 
pressed and the autocorrelation crosses zero on the scales of 
interest to us. Therefore, in order to make meaningful com- 
parisons (avoiding division by zero), the ^-series correlation 
curves were all increased by an equal, constant amount such 
that ql agrees with g6 (which has a comparable gas mass res- 
olution, but a much larger box length) at 445 km s" 1 (this ve- 
locity choice is motivated below). 

As shown in the left panels of Figure [5] the cross- 
correlation {top) and autocorrelation (bottom) from q3 agree 
well with q4 (< 3% relative difference for Av corresponding 
to less than L/4). We conclude that wvzw, which has a smaller 
gas particle mass than q3, is not significantly compromised by 
mass resolution. However, in addition to missing large scale 
power due to the L = 16 h~ l Mpc box length, the reliable sep- 
aration range (< L/4) probed by wvzw corresponds to only 
394-444 km s" 1 at z = 2-3. Conversely, the box length of 
the larger (L = 100 h~ l Mpc), but much lower-resolution (484 3 
gas particles, m gas = 9.79 x 10 7 h~ l M ^), g6 simulation should 
more than suffice. iMcDonaldl d2003l) found little difference in 
the correlation function between L = 40 and 80 h~ l Mpc sim- 
ulations. 

In Figure [6j we consider the subtracted difference between 
the correlation functions of wcw and g6 [solid lines) in order 
to characterize the effects of insufficient simulation volume 
and mass resolution and to motivate a methodology for form- 
ing hybrid correlation curves that mitigate them. Note that 
wcw is used in lieu of wvzw in order to elliminate any addi- 
tional differences due to wind models. The signature of poor 
mass resolution is illustrated (dotted lines in Fig. [6} by the 
correlation difference between ql and the mean of q3 and q4 
(which closely corresponds to the mass resolution of wcw). 
The nature of the suppression of wcw due to its small box 
length is then reflected in the residual (dash-dotted lines in 
Fig. |6| between these two curves. However, since the mass 
resolution of ql is superior to g6 by more than a factor of 2, 
the dotted line underestimates the effect for g6. Extrapola- 
tion from the right panels of Figure [5] is poorly constrained, 
although accounting for the trend implies significant flatten- 
ing of the dash-dotted line (Figure |6| on small scales. Such 
a relatively smooth alteration of the correlation function due 
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to insufficient box length is consistent with the expectation 
of constant suppressi on when the evolution o f coupled modes 
is not accounted for (McDon ald et al.l 12000, see Figure 15). 
Although the true effect is likely not a constant offset at all 
scales, this appears to be a reasonable approximation. Note 
also that the disparity between wcw and g6 at sufficiently large 
Av appears to be purely a box length effect. Therefore, we 
define the hybrid correlation function £' ! to be equal to for 
Av greater than an adopted splice velocity, v s . For Av < v„ 
£jj is equal to q vzw plus the difference between f jj 6 and q vzw 

at the splice velocity. In the case of the cross-correlation, 
this is slightly modified to preserve the boundary condition 
U (0) = £|i (0). More explicitly, 



f*(Av>v,) 
£jf(Av<v,) 

(Av<v,) 



^ 6 (Av) 

^(v s )-^(v s ) 



5 1 -ffii 



Av 



(6) 



The value of v, was chosen to be the minimum velocity 
at which the effect of the g6 mass resolution can be as- 
sumed to be negligible. In order to ensure isotropy in real- 
space, v ( must be the same for the autocorrelation and cross- 
correlation. Thus, based on Figures [5] and [6] 445 km s" 1 was 
adopted as the splice velocity (for all redshifts). It is worth 
noting that the resulting hybrid correlation curves are insensi- 
tive to the exact choice of v s due to the relative flatness of the 
difference between £ y6 and £ wvzw on these scales. 

3.2. Accounting For Spectral Resolution 

A real spectrum (i.e., observed with a telescope) is a convo- 
lution, <Sii, of the true transmitted flux along the line of sight 
with the LSF of the spectrograph. The LSF is generally Gaus- 
sian, and the width, a, determines the resolution of the data, 



j+a 



k=j-a 



2ira 



(7) 



where a must be sufficiently large with respect to a that the 
tails of the exponential are effectively zero at the limits of 
convolution. Figure [3] provides a comparison of a simulated 
spectrum at full-resolution and the same spectrum degraded 
to FWHM = 2V21n2 and a = 2.5 A. Similar to the anistropy 
introduced by redshift-space distortions, this smoothing along 
the line of sight changes the autocorrelation differently than 
the cross-correlation (Figure |7). Thus, the latter anisotropy 
must be properly accounted for in order to correct the former. 

A sensible way of determining £ CT , the correlation function 
corresponding to data of resolution a, is to smooth the simu- 
lated spectra and then compute their correlation, 
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5|| [<5„(v||.+Av,v x /),ct] 



, N, JV|| 

s «=i ii j=\ 



X 



S\\ [,5„(v||.,v ±i + Av),o-] 



(8) 



(9) 



This, however, has several disadvantages. Individually 
smoothing each spectrum is a time-consuming process which 
must be repeated for each desired value of a. Likewise, the 
correlation calculations must be duplicated, and future con- 
sideration of different resolutions requires the original sim- 
ulated spectra. More importantly, the crea tion of hybrid 
awfocorrelation curves as described in § 13.11 is only valid if 
performed at full-resolution. This is because spectral smooth- 
ing redistributes correlation along the line of sight, negating 
the validity of the splice point. Finally, smoothing the spectra 
also redistributes aliasing effects (which are mitigated in the 
hybrid correlation function) to smaller separations, limiting 
the scales which can be reliably probed at a given resolution. 
For the wvzw simulation, ±3<r (where the LSF becomes neg- 
ligible) corresponds to a third of the box length for a FWHM 
of 1.8/2.7 A at z = 2/3. 

An alternative method of accounting for spectral resolu- 
tion, applied directly to the full-resolution hybrid correlation 
function, solves each of these problems. Convolving the full- 
resolution autocorrelation function with a Gaussian LSF of 
width \/2a is mathematically identical to recalculating the 
autocorrelation with spectra smoothed by a Gaussian LSF of 
width a, 



£ll(Av),\/2<7 



= ff(Av). 



(10) 



This convenient result is due to the fact that the convolution of 
two Gaussians is itself a Gaussian and that the autocorrelation 
and spectral smoothing are both a function of radial velocity 
(see the Appendix). 

Unfortunately, this is not the case for the cross-correlation, 
and there is no corresponding analytical expression. How- 
ever, since spectral smoothing redistributes correlation along 
the line of sight, its effect in the orthogonal direction probed 
by the cross-correlation should be a relative suppression at all 
separations. The corresponding scale factor can be evaluated 
at Av = 0, where the amplitude of £J is known by virtue of 
equation[lO]and the fact that £j_ (0) = £n (0) by definition. The 
approximate solution 



ei(Av) « £l (Av) (1+0 £l( Avfl- 



where 



1 



1 



&(0) £l(0)' 



(11) 



(12) 



agrees remarkably well with results obtained using equa- 
tion [9] This is demonstrated in Figure [8] for a representative 
range of redshifts and spectral resolutions. The slight dis- 
agreement between the two methods scales with the degree of 
correlation suppression; however, the difference is < 2% for 
2 < z < 3, FWHM < 2.5 A, and B > 90". 

4. POTENTIAL SYSTEMATICS 

Simulation of the Lya flux correlation is subject to a num- 
ber of sources of systematic error. Some are either addressed 
by previous studies or may be controlled for in a limited fash- 
ion by judicious comparison of results from the simulations 
listed in Table Q] Others, we can only identify and acknowl- 
edge, but not measure or correct for. However, the primary 
interest of this study is in alterations of the correlation func- 
tion due to redshift-space distortions, for which much of this 
systematic uncertainty is mitigated. It is also worth noting 
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that while errors in the autocorrelation at small scales are pro- 
pogated to larger scales when spectral smoothing is consid- 
ered, uncertainty in the cross-correlation is only relevant at 
the scales corresponding to observed QSO pair separations 
(1.5'< 9< 4' in the case of Paper I). 

The effects of bo x len gth and mass resolution have already 
been discussed in § 13.11 For the limited scales accessible with 
the q simulation series, our hybrid correlation measurements 
appear to be largely unaffected by mass resolution and reason- 
ably well corrected for box length limitations with a constant 
offset. However, given the rapid decline of the correlation 
function on scales affected by the small box length of wvzw, 
establishing limits for the relative effect of deviations from a 
constant suppression is speculative. 

While the evolution of large scale structure at z > 2 is rela- 
tively insensitive to the cosmological parameters Vt m and VL^, 
the adopted simulation value of a% (whi ch is, however, con- 
sistent with the three-year WMAP value; ISpergel et ai]|2007l) 
likely does affect the correlation of the Lya forest. Fur- 
thermore, the simulations used in this study represent only 
a few realizations of random fluctuation amplitudes in the 
early universe. Therefore, we cannot account for any re- 
lated variance in the correlation measurements. Simulating 
the grid of amplitudes necessary for this purpose with an SPH 
code is computa tionally prohibitive at this time; however, see 
iMcDonaidl (120031) for a discussion on the alternative use of 
hydro-particle-mesh simulations. In the following subsec- 
tions, we address several remaining potential sources of sys- 
tematic error. 

4.1. Redshift Evolution 

The Az = 0.2 sampling of the wvzw simulation was used 
to verify that the redshift evolution of the autocorrelation and 
cross-correlation is smooth and well-behaved over the range 
of interest. Figure [9] illustrates this for the case of the au- 
tocorrelation (four representative Av lags are shown) at full- 
resolution with redshift-space distortions. A third order poly- 
nomial does an excellent job of fitting all seven epochs, al- 
lowing for reliable interpolation at intermediate redshifts. By 
extension, the same is assumed for the more coarsely sampled 
(Az = 0.5) g6 simulation. 

4.2. Metals 

The simulated spectra generated for this study include ab- 
sorption from H I only; however, the Lya forest in observed 
spectra is contaminated by metal lines. Associated metals can 
introduce features into the correlation function at the velocity 
difference, 

Av~cAA/A, (13) 

betwee n their absorption and th at of Lya from the same gas. 
Indeed, iMcDonald et alj (120061) found enhanced correlation 
at Av ~ 2270 km s" 1 due to Si III at rest wavelength 1206.50 
A; however, no other metal correlations were detected. More 
to the point, no metals in the IGM have known wavelengths 
closer to that of Lya than Si III; thus increased correlation 
from associated metals is not a concern for the velocity scales 
relevant to this study. Similarly, the velocity splitting of the 
Si IV doublet (Av ~ 1930 km s" 1 ) lie s bey ond our range of 
consideration, while McDonald et al. (2006) found no evi- 
dence of a correlation feature at Av ~ 500 km s" 1 correspond- 
ing to the C IV doublet. 



A third potential source of increased correlation is the clus- 
tering of metals themselves. This effect cannot be accounted 
for in the simulated spectra for two reasons. First, unasso- 
ciated metals sparsely populating the Lya forest arise from 
gas at lower redshifts, beyond the epoch for which the sim- 
ulations were run in some cases. Second, although the VZW 
wind model (see § !4.51 l has been shown to reproduce the over- 
all mass density and absorption line properties of C IV well, 
the simulations do not yet accurately reflect the clustering 
properties of metals. Fortunately, clustering of metals is not 
expected to significantly affect the Lya flux correlation func- 
tion; absorption from even the most abundant metals is 2-4 or- 
ders o f magnitude less than H I dSchave et al.l2003"tlFrve et ail 
120031 

4.3. Mean Flux Decrement Uncertainty 
The mean flux decrement (lOke & Korvcanskvll9 82). 

D A = l-(f), (14) 

where (/) is the mean of the transmitted flux (observed flux 
divided by the unabsorbed continuum flux) in the Lya forest 
can be reliably tuned to high precision in simulated spectra 
(recall § 12.3b ; however, this is only as accurate as the obser- 
vationally determined value. Measurement of Da from real 
spectra is complicated by the difficult step of estimating the 
continuum of emitted flux from the background light source 
(conveniently defined as unity in simulated spectra). For low- 
resolution spectra, the continuum has generally been extrapo- 
lated from redward of the Lya forest, assuming a power-law. 
This technique will not likely be accurate for an individual 
spectrum; however, the significant uncertainties are assumed 
to be mitigated for a sufficiently large sample. In the case 
of higher resolution spectra, a smooth continuum is fit to re- 
gions free of obvious absorption. While individually tailored, 
residual absorption will almost certainly result in artificially 
low continuum placement (corresponding to underestimated 
absorption) unless this bias can be adequately modelled. 

Numerous measurement s of Da have been made during 
the pa st two decades (see iRauchl dl99 8), Meiks in~& White! 
(2004), and references therein). The few that also determined 
its evolution as a function of redshift are compared in Fig- 
urefTUl where the thick lines represent the redshift range of the 
data u sed. Consistent with the above discussion, [Press et all 
( 1 1993b extrapolated the continuum for 29 quasars and ob- 
tained 

Df 3 (z)=l- e - 00037 (' +z » 346 , (15) 

whereas continu um fits to echelle - resolu tion spectra by 
iKim et alJ (1200 ll) and iKirkman et all (120051) yielded signifi- 
cantly lower values (less absorption). The latter found 

£»f 05 (z) = 0.0062 (1+z) 2 - 75 (16) 

and claimed errors of less than 1 % b ased on test s using ar- 
tificial spectra. A reevaluation of the iPress et aTl (1 19931) re- 
sults bv lMeiksin & W hite (2004) reported mu ch better agree- 
ment w ith the high resolution studies; however. lBernardi et alj 
(120031) similarly extrapolated the continua for a sample of 
1061 quasars and produced results very similar to the original 
IPress et ail d 19931) values. The mean flux decrement remains 
observationally uncertain. 

Unfortunately, as has been shown previously for mat- 
ter power spectrum m e asurements made with t he Lya for- 
est dCroft et alj l2002bt IZaldarriaga etaTJ 120031: ISeliak et alj 
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2003), both the amplitude and shape of the correlation func- 
tion are sensitive to Da. Figure [TT] shows the percent differ- 
ence in the correlation function for simulated spectra tuned 
to have the mean flux decrement prescribed by either equa- 
tion[l5]or equation [TBI In the bottom panel of FigureQT] the 
correlation functions compared have been arbitrarily scaled to 
unity at 400 km s" 1 in order to mitigate differences solely in 
amplitude. We have addressed this systematic uncertainty by 
carrying o ut our analy s is usin g th e mean flux decreme nt val- 
ues of both lPress etal] (1 19931) and lKirkman et all d2005l) . Un- 
less otherwise stated, results from the former are used in the 
figures throughout this paper (where this choice is secondary 
to other effects being considered). 

4.4. Spectral Resolution 

Section [3~2l discussed how to account for arbitrary spectral 
resolution when using the correlation functions computed at 
full-resolution. Here we consider in greater detail the sen- 
sitivity of the correlation function to small changes in spec- 
tral resolution (or uncertainty in that parameter). Figure Q~2] 
shows the relative difference in autocorrelation correspond- 
ing to a 4% change (AFWHM = 0. 1 A for FWHM = 2.5A) in 
resolution. This difference scales roughly linearly for larger 
AFWHM and is less for the cross-correlation. Thus, treating 
data with FWHM = 2.3 and 2.5 A as having the same resolu- 
tion introduces an error of up to 5% - 6%. 

4.5. Wind Model 

Prescriptions for galactic winds, which transport pro- 
cessed gas from within galaxies to the surrounding IGM, 
are relatively new a d dition s to cosmological simulations. 
Sprin geT& HernquisH (120031) incorporated a constant wind 
(cw) model in order to reduce the amount of gas available for 
star formation in galaxies. Essentially, a fraction of the gas 
particles, dictated by the current star formation rate and a rel- 
ative mass loading factor r\, are ejected from a galaxy via su- 
perwinds. They then travel without hydrodynamic interaction 
at a constant velocity v lr m</ until the SPH density falls below 
10% of the critical density for multi-phase collaps e. Based 
on earlier simu lation work b y lAguirre et alj (1200 1|) and ob- 
servations from lMartinl(ll999l) and lHeckman et al.1 d2000l) . the 
two free model parameters were set at v,„, 1£ / = 484 km s" 1 and 
r\ = 2. This yields broad agreement with observations of the 
stellar mass density at z=0; however, th e wind velocity is un- 
physic ally large for small galaxies, and lOppenheimer & Pavel 
(2006) found that C IV is overproduced in the IGM com- 
pared to observed Slciv data. These authors also note that the 
CW model does not converge well with resolution. That is, 
in higher-resolution simulations which resolve small galaxies 
earlier, the winds turn on earlier and heat the IGM in excess 
of lower-resolution simu la tions. 

Opp enheimer & Pavel (120061) also investigated several, 
more sophisticated prescriptions for galactic outflow, con- 
trasting their effect on the IGM and comparing the results to 
observational data. The most successful models were vari- 
ants of momentum-driven winds. In the case of VZW, the 
wind speed, v w m = 3 cry r fj~— 1, and the mass loading fac- 
tor, r\ = cr„cr~ l , both scale as the galaxy velocity dispersion, 
g = -J— $ /2. Here, $ is the gravitational potential, and 
f L = f LQ x io 0029 ('«« z+9 ) 25+0 - 417694 is the galaxy luminosity 
in units of its critical luminosity. The free parameter a a was 
chosen to be 300 km s" 1 , corresponding to a Salpeter initial 
mass function and a typical starburst spectral energy distri- 



bution, and //, i0 was al lowed to vary rando mly in the range 
1.05-2 as observed by iRupke etal] (120051) . Unlike the cw 
wind model, VZW was shown to non-trivially reproduce a 
wide range of C IV absorption observations. 

While more detailed studies of the effec ts of galactic winds 
on the Lyq forest have been carried out (ICroft et al.l l2002a: 
iPesiacques et al.l 120041 iMcPonald et al.ll2005l) . our primary 
interest is in investigating how the different wind models in- 
cluded in the g6 (cw) and wvzw (VZW) simulations might af- 
fect our Lya forest flux correlation measurements. The wcw 
and wnw simulations are identical to wvzw with the excep- 
tion of their wind models. As their nomenclature indicates, 
the former incorporates the cw model, while the latter in- 
cludes no winds at all (NW). Slight differences in the flux 
distribution caused by the inclusion of winds were mitigated 
by the rescaling of opacities described in § 12.31 Although the 
three simulations are identically affected by box length limi- 
tations, corrections were applied as described in § 13. ll so that 
meaningful comparisons could be made of correlation curves 
that otherwise cross zero in the region of interest. Figure Q~3] 
shows the percent difference between the correlation values 
obtained from each of these two simulations and those from 
wvzw (in redshift-space, at z = 3, and at full-resolution). The 
wcw and wvzw results differ by < 1 % and < 4%, for the cross- 
correlation and autocorrelation respectively, on scales larger 
than ~ 50 km s" 1 , indicating that our correlation measure- 
ments for g6 and wvzw are only marginally affected by the use 
of different wind models. Furthermore, while we assume that 
inclusion of the currently preferred wind model yields more 
accurate results than neglecting galactic winds altogether, the 
wnw and wvzw comparison demonstrates that these two ex- 
tremes represent a difference of only < 7%. 

5. IMPLICATIONS FOR THE AP TEST 

5.1. Signal-To-Noise Ratio 

Cross-correlation measurements were repeated with vary- 
ing degrees of Gaussian noise added to the individual simu- 
lated (wvzw) spectra. Although this has no effect on the mean 
correlation values (which have been averaged over many lines 
of sight), signal-to-noise ratio (S/N) does affect the disper- 
sion of those values. However, even for relatively low S/N, 
the corresponding increase in is negligible relative to the 
intrinsic variation in £ between different lines of sight. The 
latter scales inversely with path length, but even for the entire 
Lya forest redward of Ly/3 absorption, the difference in 
between S/N = 5 and S/N = oo is < 2% for 2 < z < 3, resolu- 
tion FWHM < 2.5A, and < 9 < 300 arcseconds This is in 
agreement with the assertion by iMcPonaldl (120031) that only 
moderate quality data is needed for a large number of quasar 
pairs to carry out the Alcock-Paczynhski test. The S/N re- 
quirements of observed spectra are dictated not by correlation 
measurements, but by the needs of reliable continuum fitting. 

5.2. Continuum Errors and Da Variance 

Errors in fitting the continua of observed QSO spectr a can 
affect calculation of the mean flux decrement (recall § I4.31 l 
as well as correlation measurements. Comparison of Da for 
a particular spectrum to the expected mean flux decrement 
might be used, in principle, to constrain systematic errors in 
the determination of the continuum. However, genuine vari- 
ation in Da arises naturally between lines of sight (decreas- 
ing with increasing path length) due to finite sampling of the 
local large scale structure. Simulated spectra provide an op- 
portunity to quantify the expected distribution of mean flux 
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decrement measurements in the absence of continuum fitting 
errors. 

Table [2] provides the variance in Da as a function of red- 
shift and path length (in units of h~ l comoving Mpc) for wvzw 
and g6. Although some validation is given by the general 
agreement between the two simulations, the g6 results are sys- 
tematically lower than those for wvzw (the percent difference 
increases from approximately 1% to 9% at z = 3 and 2, re- 
spectively). If the difference was dominated by the diversity 
in large scale structure contained within the different simula- 
tion volumes, one would expect the g6 variances to be larger. 
Since this is not the case, we presume that the differences pri- 
marily reflect the greater mass resolution of the wvzw simula- 
tion (note that this is consistent with the difference increasing 
monotonically as the fraction of pixels in low density regions 
increases at lower redshift). At z = 2.2, the standard deviation 
in Da for a path length of Az = 0.2 is <jjy A — 0.017, correspond- 
ing t o 9.0% and 12 1 % of th e D A value from lPress etail (119931) 
and iKirkman et alJ (120051) . respectively. This decreases to 
<j Da ~ 0.010 (5.5% and 7.3%) for Az = 0.545, the path length 
of the full "pure" Lya forest (redward of the onset of Ly(3 
absorption). 

5.3. Anisotropy Corrections 

The primary goal of this work is to model anisotropies in 
the observed Lya forest correlation function, facilitating a 
new application of the AP test using spectra of QSO pairs 
(such as those presented in Paper I). To this end, we have com- 
puted the autocorrelation and cross-correlation in both real 
and redshift-space, investigated potential sources of system- 
atic error, and considered the impact of spectral smoothing. 
Our full-resolution, hybrid correlation measurements are pro- 
vided in Tables [3]|6] (complete versions of the stubs included 
here can be found in the electronic edition of ApJ or upon re- 
que st) for the mean flux de crements of bofh lPress et alfd 19931) 
and IKirkman et al] (2005). Note that the velocity scales are 
redshift dependent, so a unitless parameterization (first col- 
umn) is used which is not the same for the autocorrelation 
and cross-correlation. 

Implementation of the AP test itself is nontrivial and the 
subject of Paper III in this series. However, we conclude by 
outlining a scheme for the use of these simulation results that 
mitigates the systematic uncertainty discussed in § |4j To re- 
iterate, the correlation function of the resolved Lya forest is 
isotropic in real-space, and adjusting the angular diameter dis- 
tance until cross-correlation measurements (the data) agree 
with the autocorrelation (the model) yields the correct cos- 
mology. 

Figure [14] shows the effects of redshift-space distortions 
and spectral smoothing on the cross-correlation (top left). 
These can be accounted for in observed cross-correlation 
measurements by applying the ratio of the full-resolution, 
real-space simulated cross-correlation divided by its counter- 
part for smoothed data (using equations[TTland[T2li in redshift- 
space (bottom left panel of Figure fT4l) . This correction re- 
quires adopting an angular diameter distance and, therefore, 
must be applied independently for each cosmology consid- 
ered. Using the ratio of simulation results allows for partial 
cancellation of systematic errors. The right panels of Fig- 
ure Q3] show the relative difference in th ese correcti o ns be- 
tw een using the mean fl ux decrement of iPress et al.l (1 19931) 
or IKirkman et al.1 (120051) . While still a significant source of 
systematic uncertainty, the impact of Da on the correlation ra- 
tio is reduced relative to the correlation function itself (recall 



Figure [TTJ. 

Until sufficient high-resolution (echelle) data exists for reli- 
able determination of the autocorrelation (many lines of sight 
are needed to compensate for significant variance), simulated 
measurements provide the only reasonably continuous model. 
However, more abundant observational data obtained at lower 
resolution can be used to correct systematic error in the sim- 
ulation data. This is accomplished by smoothing the full- 
resolution, redshift-space correlation curve (recall that the 
Z = 3 hybrid autocorrelation is shown in Figures |4] and [JJ as 
appropriate (using equations 171 andfTUli and fitting it to the ob- 
served data. The same corrections can then be applied to the 
simulated full-resolution, real-space autocorrelation model, 
which is not affected by the discussed anisotropies. 

6. SUMMARY 

Using cosmological hydrodynamic simulations, we have 
modelled the Lya flux autocorrelation and zero-lag cross- 
correlation in both real-space and redshift space at 1.8 < 
z < 3. Mock Lya flux absorption spectra were generated 
from eight SPH simulations with and without inclusion of 
redshift-space distortions caused by Hubble expansion, bulk 
flows, and thermal broadening. The simulations considered 
(wl6n256vzw at 1.8 < z < 3.0, g6 at 1.5 < z < 3.0, and 
wl6n256cw, wl6n256nw, ql, q2, q3, and q4 at z = 3) primar- 
ily differ in their size, mass resolution, and prescription for 
galactic outflow. The lines of sight through each simulation 
box were selected such that different pairings form 73 unique 
transverse separations spanning the range 0-5 arcminutes. 
Our analysis is summarized below. 

1) Autocorrelation and zero-lag cross-correlation measure- 
ments were computed from the extracted spectra for both real- 
space and redshift-spa ce and for the mean flux decrement val- 
ues re ported by both IPress et all ( 1 1993b and IKirkman et al.l 
(2005). The difference in the autocorrelation and cross- 
correlation corresponding to this observationally uncertain 
parameter was found to be 20% -45% and 20% -35%, re- 
spectively, affecting both the shape and amplitude. 

2) Convergence of the simulated Lya flux correlation as 
a function of mass resolution was tested at z = 3 for the 
Gadget code using the ^-series simulations (which identi- 
cally evolve different numbers of particles within boxes of 
equal volume). The difference in autocorrelation and cross- 
correlation between q3 (m gas = 3.72 x 10 6 h~ l M Q ) and q4 
(m gas = 1.10 x 10 6 h' 1 Mq) is less than 3% on all scales. 

3) The g-series was also used to characterize the effect of 
insufficient mass resolution in g6 and, indirectly, the effect 
of the inadequate simulation volume of wl6n256vzw. In or- 
der to correct for these limitations of current simulations, hy- 
brid correlation curves were then formed by splicing together 
those from wl6n256vzw and g6 at Av = 445 km s . At 
smaller velocities, the hybrid correlation is equal to that of 
wl6n256vzw plus a constant boxsize correction (in the case 
of the cross-correlation, this is slightly modified to preserve 
the boundary condition at Av = 0). At larger velocities, where 
the effects of mass resolution were projected to be insignifi- 
cant, the hybrid correlation is provided by that of g6 without 
alteration. 

4) An approximate solution is presented for obtaining the 
zero-lag cross-correlation corresponding to arbitrary spectral 
resolution directly from the zero-lag cross-correlation com- 
puted at full-resolution (an exact solution is available in the 
case of the autocorrelation). This approximation is good to 
within 2% for the relevant redshift range at velocity differ- 
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ences corresponding to angular separations greater than 90 
arcseconds. 

5) The effects of three prescriptions for galactic out- 
flow on the Lya flux correlation were investigated with the 
wl6n256vzw, wl6n256cw, and wl6n256nw simulations. The 
difference between the preferred variable-momentum wind 
model (vzw, used for wl6n256vzw) and the older constant 
wind model (CW; used for g6) was found to be < 1% and < 
4% at scales larger than « 50 km s" 1 for the cross-correlation 
and autocorrelation respectively. The corresponding differ- 
ence between VZW and no winds at all increases to only < 5% 
and < 7%. 

6) For an adopted mean flux decrement, the variance from 
one line of sight to another was computed as a function of 
redshift and path length. At z = 2.2, the standard deviation in 
Da for a path length of Az = 0.2 is <jd a — 0. 017, corresponding 
to 9 0% and 12.1% of th e D A value from iPress etail d 19931) 
and Kirkm an" et alj (|2005), respectively. 

7) Aside from those sources of systematic error already 
summarized above, we find that redshift evolution of the Lya 
flux correlation is sufficiently sampled for reliable interpola- 
tion and argue that absorption from metals is insignificant. 
The evolution of large scale structure at z > 2 is not sensi- 
tive to the values for the cosmological parameters Q m = 0.3 
or fl\ = 0.7 assumed by the simulations considered here, and 
as = 0.9 is consistent with the three-year WMAP value. Sys- 
tematic error associated with variance of random fluctuation 



amplitudes in the early universe or deviations from a constant 
offset due to finite boxsize cannot be addressed with currently 
available simulations. 

8) Correcting for anisotropies due to redshift-space dis- 
tortions and spectral smoothing with ratios of the correla- 
tion measurements allows for significant reduction in system- 
atic error. The maximum difference betwe en using the mean 
flux de crements of either Pres s et all d 1993b or lKirkman et alj 
(2005) (the dominant source of uncertainty) decreases to 
8%- 16% at 2 < z < 3, and presumably the true value is in- 
termediate. We describe a simple scheme for implementing 
our results, while mitigating systematic errors, in the context 
of a future application of the AP test using observations of the 
Lya forest in pairs of QSOs. 
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FIG. 3. — Simulated spectra from wvzw and g6 show the larger box length (path length) of the latter, the higher mass resolution (narrower features) of the 
former, the subtle redistribution of opacity due to redshift- space distortions, and the smoothing effect of spectral resolution. 




FIG. 4. — In real-space, the correlation function is isotropic. However, in redshift-space, distortions caused by line-of-sight velocities affect the autocorrelation 
(£||) and cross-correlation (£x) differently. 
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FIG. 1 1 . — Amplitude and shape of the Ly« fl ux correlatio n function are s ensitive to the mean fl ux decrement, as evidenced by the relative difference in £ for 
simulated spectra tuned to match the values from Press et al. 1 1993, P93) and Kirkman et al. (2005, K05). 
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Relative difference in autocorrelation corresponding to a 4% change in spectral resolution FWHM. 
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FIG. 13. — Percent difference in correlation between wvzw and wI6n256[cw/mv] indicates that Lya correlation measurements are marginally affected by the 
prescription for galactic outflow. The preferred momentum-driven wind model VZW differs from the older constant wind model CW at z = 3 and on scales larger 
than ~50 km s" 1 by < 1% for the cross-correlation and < 4% for the autocorrelation. This increases to < 5% and < 7%, respectively, when compared to no 
winds at all (NW). 
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TABLE 1 
Simulation Properties 



Simulation 




N p b 






Wind 


Redshift 




Ali 


Name 


Alias 


(A -1 Mpc) 




(r'M ) 


(/r 1 kpc) 




Range 


Az 




ql 




10 


2x64 3 


4.24 x 10 7 


6.25 


cw 


3.0 




30000 


q2 




10 


2x96 3 


1.25 x 10 7 


4.17 


cvv 


3.0 




30000 


q3 




10 


2 X 144 3 


3.72 x 10 6 


2.78 


cw 


3.0 




30000 


q4 




10 


2X216 3 


1.10 x 10 6 


1.85 


cvv 


3.0 




30000 


wl6n256nw 


wnw 


16 


2 x 256 3 


2.71 x 10 6 


1.25 


NW 


3.0 




20000 


wl6n256cw 


wcw 


16 


2 x 256 3 


2.71 x 10 6 


1.25 


cvv 


3.0 




20000 


wl6n256vzw 


wvzw 


16 


2 x 256 3 


2.71 x 10 6 


1.25 


vzw 


1.8-3.0 


0.2 


20000 


g6 


100 


2 x 484 3 


9.79 x 10 7 


5.33 


cvv 


1.5-3.0 


0.5 


3000 


Note. — 


The columns (left to 


right) are: 


simulation name and abbreviation, box length, total number of particles, s 


;as particle mass, 


equivalent Plummer 


gravitational softening length, wind model, redshift range and interval, and number of sets of lines of sij 


jht extracted alonj 


; each of the three principal 


axes. 



a These physical scales are given in comoving coordinates. 
b The total number of particles is evenly divided between dark matter and gas. Therefore, the dark matter particle mass, m c i m , is simply the gas particle mass 
scaled by (f! m — f2j)/Qj. 

c The simulations have either no prescription for galactic outflow (NW), a constant wind (CW) model, or momentum-driven winds (VZW). See jj |4.5l for details. 



TABLE 2 

D A Variance 



z 








wvzw 






1.5 






0.024144 


1.8 


0.041109 






2.0 


0.048716 




0.044152 


2.2 


0.056126 






2.4 


0.063490 






2.5 






0.064115 


2.6 


0.069438 






2.8 


0.074085 






3.0 


0.076494 




0.075647 



a The variable / is the path length in h 1 comoving Mpc 



TABL E 3 

Cross-Correlation JPress et al.I dT993h Mean Flux Decrement) 



Av a 






10x£j 


(Av) in real-space 










10x£j 


(Av) in z-space 






(<5v x (z)) 


2=1.8 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


z = 1.8 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


0000 
0001 
0002 
0003 

I 

0300 


0.4315 
0.4263 
0.4151 
0.3997 


0.5914 
0.5853 
0.5702 
0.5493 


0.7994 
0.7917 
0.7718 
0.7436 


1.0641 
1.0541 
1.0275 
0.9903 


1.4036 
1.3907 
1.3555 
1.3067 


1.8199 
1.8026 
1.7566 
1.6933 


2.3333 
2.3094 
2.2504 
2.1686 


0.7489 
0.7302 
0.7089 
0.6870 


0.9812 
0.9614 
0.9364 
0.9090 


1.2594 
1.2392 
1.2103 
1.1772 


1.5876 
1.5660 
1.5316 
1.4915 


1.9748 
1.9499 
1.9089 
1.8593 


2.4281 
2.3995 
2.3508 
2.2893 


2.9558 
2.9216 
2.8630 
2.7876 


0.0078 


0.0107 


0.0176 


0.0230 


0.0238 


0.0237 


0.0461 


0.0489 


0.0555 


0.0654 


0.0840 


0.1030 


0.1158 


0.1270 


Note. — TIickl- 


arc ihe hybrid correlation values 


for bill-resolution. [The complet 


version of this lable can be found 


n ihe electronic 


edition of ApJ or upon request.] 













a 74 Av values given in units of <5vj_ (z) = 0. 67455; + 0. 38896 km s ' (or, alternatively, arcseconds) 
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TABLE 4 

Autocorrelation <Press et AL.N1993D Mean Flux Decrement) 



(<5v||(z)) 






10 X £|| (Av) in real-space 










10 


II (Av) in z-space 






2=1.8 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


z=1.8 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


0000 


0.4315 


0.5914 


0.7994 


1.0641 


1.4036 


1.8199 


2.3333 


0.7489 


0.9812 


1.2594 


1.5876 


1.9748 


2.4281 


2.9558 


0001 


0.4250 


0.5837 


0.7902 


1.0530 


1.3902 


1.8036 


2.3135 


0.7478 


0.9797 


1.2572 


1.5847 


1.9709 


2.4231 


2.9493 


0002 


0.4097 


0.5649 


0.7672 


1.0249 


1.3558 


1.7616 


2.2622 


0.7445 


0.9750 


1.2507 


1.5759 


1.9593 


2.4081 


2.9303 


0003 


0.3897 


0.5398 


0.7359 


0.9862 


1.3079 


1.7027 


2.1899 


0.7392 


0.9674 


1.2400 


1.5615 


1.9404 


2.3838 


2.8994 


I 

0999 


7.5E-4 


6.6E-5 


-0.0019 


-0.0023 


0.0014 


-0.0098 


0.0026 


-0.0028 


-0.0042 


-0.0060 


-0.0078 


-0.0094 


-0.0116 


-0.0114 


Note. — These are the hybrid t 




s for full-resoluii 


on. [The eomplel 


; version of this table can be found 


in the electronic 


; edition of ApJ or upon request.] 













3 1000 A vvalues given in unils of 5vi| (z) = 0.25246 z+ 1 .473] kms 



TABLE 5 

Cross-Correlation (Ki rkmanet al.I ( 120051) Mean Flux Decrement) 



Av a 






lOx^j 


( Av) in real-space 










10x£j 


(Av) in z-space 






(Sv±(z)) 


z=1.8 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


z = 1.8 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


0000 
0001 
0002 
0003 

I 

0300 


0.3219 
0.3173 
0.3078 
0.2950 


0.4503 
0.4449 
0.4319 
0.4143 


0.5798 
0.5735 
0.5569 
0.5340 


0.7403 
0.7323 
0.7110 
0.6817 


0.9426 
0.9326 
0.9052 
0.8678 


1.1868 
1.1738 
1.1389 
1.0918 


1.4890 
1.4715 
1.4279 
1.3684 


0.5778 
0.5596 
0.5400 
0.5204 


0.7721 
0.7527 
0.7295 
0.7049 


0.9534 
0.9342 
0.9081 
0.8793 


1.1649 
1.1448 
1.1144 
1.0807 


1.4131 
1.3905 
1.3555 
1.3151 


1.7022 
1.6772 
1.6371 
1.5886 


2.0396 
2.0107 
1.9642 
1.9065 


0.0057 


0.0074 


0.0119 


0.0148 


0.0139 


0.0119 


0.0274 


0.0380 


0.0400 


0.0443 


0.0570 


0.0707 


0.0786 


0.0870 


NOTE. — These 


are the hybrid ce 


(relation values 


for full-resolulion. [The eompleli 


: version of this tabli 


: can be found i 


n the electronic edition of ApJ or upon request.] 













3 74 Av values given in unils of 5vj_ (z) = 0. 67455; + 0. 38896 km s * (or, alternatively, arcseconds) 



TABLE 6 

Autocorrelation ( Kirk man et al.I j2005l) Mean Flux Decrement) 



Av a 






10 X i; || (Av) in real-space 










10 x<; 


II (Av) in z-space 






(5v||(2» 


2=1.8 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


2=1.8 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


0000 


0.3219 


0.4503 


0.5798 


0.7403 


0.9426 


1.1868 


1.4890 


0.5778 


0.7721 


0.9534 


1.1649 


1.4131 


1.7022 


2.0396 


0001 


0.3162 


0.4435 


0.5721 


0.7312 


0.9320 


1.1743 


1.4742 


0.5770 


0.7709 


0.9518 


1.1628 


1.4103 


1.6987 


2.0352 


0002 


0.3033 


0.4275 


0.5531 


0.7087 


0.9052 


1.1425 


1.4362 


0.5746 


0.7674 


0.9471 


1.1565 


1.4022 


1.6883 


2.0221 


0003 


0.2868 


0.4065 


0.5277 


0.6782 


0.8685 


1.0985 


1.3834 


0.5707 


0.7616 


0.9393 


1.1463 


1.3890 


1.6715 


2.0007 


I 

0999 


6.4E-4 


3.2E-5 


-0.0011 


-0.0015 


0.0011 


-0.0047 


0.0012 


-0.0021 


-0.0030 


-0.0041 


-0.0054 


-0.0067 


-0.0083 


-0.0084 


Note. — These 


are the hybrid correlation vaLics for full -resolution. [The complel 


aversion of thi: 


» table can be found 


in the electronic edition of ApJ or i 


ipon request.] 













a 1000 Av values given in unils of 5vn (z) = 0.25246;+ 1.4731 km s" 
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APPENDIX 

AUTOCORRELATION CALCULATION FOR ARBITRARY SPECTRAL RESOLUTION 

One way to compute the simulated correlation function corresponding to a given spectral resolution is to carry out the calcula- 
tions using spectra which have been individually smoothed as appropriate. However, assuming a Gaussian line spread function 
(LSF), the autocorrelation curve for data of arbitrary spectral resolution can also be obtained by simply convolving the full reso- 
lution curve with the LSF broadened by a factor of a/2 (eqs.|H|7][8] andfTOb. The validity of this relation for our discrete, periodic 
simulated spectra has been tested and verified. Here, in the interest of clarity, we demonstrate its origin for the simplified case of 
continuous spectra. In this limit, equations [4][7] and[8]become 



i 11 r 

£n ( Av) = jf E E t^t / s " ( v h v ±) s » ( v \\ + Av > v ±) dv . 



5n 



f(r) 



1 ±irl 

e 2a 2 dr, 



2ir<j 



and 



, N N ± 
s ll v ' V. /L^i M, ^—^ 



respectively. Combining (IA21 > and (lA3t . and rearranging, yields 

n=\ 

After substituting = v 



Sn (v'|| , Vj_ ; ) , cr 



i A i N± 



S ll V ' M IV, £j 



S„(j,V_Li) S„(f,v ±i ) 



S n (v||+Av,v_l,) ,<j 



dv\\ 



7.-KG 1 



■e e 



dv\\ drdf. 



t and noting that ^-^e^e ^ 2 dp= „ j= e *^ , this becomes 

1 



N± 



11= 1 (=1 J 1 



5n(T,V±i) S n (f,Vj_i) 



2y / 7T<7 



e i" 1 drdf. 



With another change of variables (f = r + a), £° can now be written in terms of ( IAU and dA2b 

1 



1 1 N± 1 

t(Av)=^E^E ra 

»=1 i=l J ii 



ef|(Av) 
|[f(Av)=5| 



' 1 1 N± 1 

a^E^E jj V|| 



2y / 7T(7 

4 (t, v_l,) 4 (r + a, v_u) rfr 



-(Q-Av) Z 

e drda 



1 



27t(V2ct) 



e 2(V2t) 2 qPq; 



|||(Av),\/2cr 



(Al) 
(A2) 

(A3) 
(A4) 



(A5) 



(A6) 

(A7) 
(A8) 



